Mass, Charge and Motion in Covariant Gravity Theories 
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Previous work established a universal form for the equation of motion of small bodies in theories 
of a metric and other tensor fields that have second-order field equations following from a covariant 
Lagrangian in four spacetime dimensions. Differences in the motion of the "same" body in two 
different theories are entirely accounted for by differences in the body's effective mass and charges 
in those different theories. Previously the process of computing the mass and charges for a particular 
body was left implicit, to be determined in each particular theory as the need arises. I now obtain 
explicit expressions for the mass and charges of a body as surface integrals of the fields it generates, 
where the integrand is constructed from the symplectic current for the theory. This allows the entire 
prescription for computing the motion of a small body to be written down in a few lines, in a manner 
universal across bodies and theories. For simplicity I restrict to scalar and vector fields (in addition 
to the metric), but there is no obstacle to treating higher-rank tensor fields. I explicitly apply the 
prescription to work out specific equations for various body types in Einstein gravity, generalized 
Brans-Dicke theory (in both Jordan and Einstein frames), Einstein-Maxwell theory and the Will- 
Nordvedt vector-tensor theory. In the scalar-tensor case, this clarifies the origin and meaning of the 
"sensitivities" defined by Eardley and others, and provides explicit formulae for their evaluation. 



While general relativity (together with a cosmological 
constant) appears to provide a successful macroscopic de- 
scription of all known gravitational phenomena, it is of 
interest to explore alternative theories that may provide 
a more fundamentally appealing description or suggest 
new experiments leading to the discovery of new phe- 
nomena. A basic framework for inventing and analyzing 
such theories is provided by Lagrangian field theory of a 
Lorentz metric and other tensor fields, and indeed some 
version of this framework is applied in essentially all mod- 
ern attempts. In the present work we will employ this 
framework and restrict further to theories that are dif- 
feomorphism covariant (no background structure), have 
second-order field equations, and live in four spacetime 
dimensions. This class contains a majority of the theories 
commonly considered, and many theories that apparently 
violate these criteria may be rewritten so as to satisfy 
them. For example, it is well known that f(R) theories 
(with fourth-order field equations) may be rewritten as 
scalar-tensor theories with second-order field equations, 
and many higher-dimensional theories are analyzed with 
a four-dimensional effective Lagrangian. 

Previous work [l[ (hereafter paper I) considered the 
problem of the motion of small bodies in such theories. 
Here, a "small body" is defined in terms of its external 
fields: if a region surrounding the putative body may 
be found such that each field approximately splits into a 
1/r-type "body field" plus a smoothly varying "external 
field" , then one says that a small body is indeed present. 
A perturbative formalism makes this notion precise, and 
a zero-size limit identifies the body region with a world- 
line in a smooth background spacetime. We found in 
paper I that the equation for this worldline takes a form 
completely independent of both the details of the body 
composition and the details of the Lagrangian governing 
the theory. The only parameters that appear are cer- 



tain effective mass and charges of the body, which may 
be determined for a particular body in a particular the- 
ory by a process that involves the solution of the field 
equations for that body with boundary conditions repre- 
senting the state of the nearby "external universe" . In 
paper I, the details of this computational process were 
left implicit, to be determined in each specific theory as 
the need arises. The main new result of the present work 
is an explicit formula for the effective mass and charges 
that holds independent of the body and theory. 

The formula is given in an abstract form in equation 
(|88[). below, and is expressed for a scalar- vector-tensor 
theory in equations (|120[) - (ll23p . The main ingredient 
is the symplectic current uj a , which follows from a sec- 
ond variation of the Lagrangian, equations (JTUJ) and . 
The symplectic current depends on two variations (per- 
turbations) of the fields, and thus may be thought of as 
a function with two "slots" . Our formulae instruct one 
to plug in to these slots 1) the 1/r fields generated by 
the body (denoted a/r) and 2) a certain "test variation" 
/ that controls which charge one computes. The 1/r 
fields are to be computed in the near-zone, i.e., from the 
far-field of a solution of the field equations for an iso- 
lated body. One then plugs this portion of the field in 
to the formulae to determine the charges of that body. 
In their general form these formulae involve a surface in- 
tegral taken over a vanishingly small sphere centered on 
the position of the body in the background spacetime. 
However, in any specific theory it is possible to rewrite 
the integral in near- zone form, i.e., over a large sphere 
surrounding a finite-size body. This bypasses the need to 
explicitly single out the 1 Jr terms and provides a surface 
integral that may be applied directly on the near-zone 
fields. We provide such a formula for each specific theory 
we consider. 

A simple example that indicates the nature of our 
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results is the classic Brans-Dicke theory analyzed 
in the Jordan frame. In all scalar-tensor theories of 
our class, the motion of a body through a background 
{gab, 4 1 } is given in terms of (non-constant) parameters 
M (mass) and q (scalar charge) by 



Mu b \7 b u a = q (g ab + u a u b ) V b 



(1) 



where u a is the four- velocity and everything is evaluated 
on the worldline. To apply this formula in Brans-Dicke 
theory (or any specific theory) , one must determine q and 
M for the body of interest. Our formalism instructs one 
to take a near-zone viewpoint, where the body is of finite 
size and immersed in an external universe whose scale of 
variation is large compared to the body. In particular, 
denote the value of <p on the worldline at the present time 
by 4>. Then, one should consider stationary solutions of 
the field equations with asymptotic behavior — > f]^ 
and (p — > <fi, and containing a body of the type under 
consideration. Having obtained such a solution, the as- 
sociated charges are given by 

M = -L [ r 2 dn n l [cp (d k g kl - d l9kk ) - 28^} (2) 
16tt J x 



-1 
16tt 



q = - — / r 2 dfl n 



2-di<j) + dig 00 + d k g ki - dig kk 



(3) 



where r = ^8ijX t xi , n l = x l /r, dVL is the area ele- 
ment on the unit sphere, repeated indices are summed, 
and J indicates an r — > oo limit of the integral^ 
Given a stationary, asymptotically constant solution of 
the field equations, equations © and ([3]) provide surface- 
integral formulae for the mass and charge in the spirit 
of the Arnowitt-Deser-Misner (ADM) formula for the 
mass in general relativity pj. In particular, they are 
coordinate-invariant under the restriction of asymptot- 
ically Minkowskian metric components. However, our 
mass M is not derived from Hamiltonian considerations 
or related to time-translation symmetry in any way; in- 
stead, we simply refer to M and q as the "mass" and 
"charge" given their appearance in equation for the mo- 
tion of the body, equation ([TJ . 

To provide a deterministic prediction for the motion of 
a body, it is necessary not only to determine the charges 
M and q at an initial time when the body configuration 
is prescribed, but to have a prescription for determining 
M and q at any given time. The simplest way of doing 
so is to choose enough properties of the body (such as its 
shape and total baryon number) such that the near-zone 
solution becomes a deterministic function of the asymp- 
totic scalar field value <f>. In this case the mass and charge 



likewise become functions of <j> by equations ([2]) and ([3]), 
and these functions may be plugged into equation ([TJ to 
determine the motion. It is important to note, however, 
that not all prescriptions one might try will be consistent 
with the theory. In particular, we find, with the same sta- 
tus as equation (JXJ) (i.e., holding in a general theory), an 
equation for the evolution of the mass, 



u b V b M 



-qu b V b <t>. 



(4) 



Under circumstances where the body is assumed a deter- 
ministic function of the external field value, equation Q 
implies that the functions M (0) and q(<fi) must satisfy 
M'((f) = —q(<f>). This relates our scalar charge q to the 
sensitivity of the mass to changes in the external scalar 
field value, reproducing a basic feature of the classic ap- 
proach of Eardley @ . (Indeed, our results are completely 
consistent with previous work; see discussion in section 
[Vl below) . We go beyond previous work by treating gen- 
eral theories in a unified framework, by increasing the 
rigor of the derivation, by allowing bodies whose field 
configuration is not a deterministic function of the exter- 
nal field value, and by providing explicit surface-integral 
formulae for the evaluation of the mass and charge(s). 

In some references the equation q — —M'(<j>) is viewed 
as defining the scalar charge q (or equivalently the sen- 
sitivity s = —q<j)/M). This definition is less general than 
ours, since it requires a family of solutions (parameter- 
ized by <j>), whereas our q can be computed for a single 
solution from equation ([3J (or from analogous equations 
in more general theories). However, the definition is also 
incomplete in the sense that no instructions are given as 
to what property of the body is to be held fixed when 
the derivative with respect to $ is takend Our viewpoint 
turns this question around, stating that one should find 
a stationary solution, compute q from equation ([3]) and 
use the relationship q = —M'(cj>) to help determine what 
properties of the body are held fixed under changes in cf>. 
In general, the only foolproof method of determining how 
a body responds to such changes is to solve the field equa- 
tions for that body immersed in a slowly time-dependent 
external field. However, our equation indicates a par- 
ticular property that any such solution will satisfy: the 
mass M, as defined by equation |(3J), will change accord- 
ing to the scalar charge q, as defined by equation ([3j, by 
M'(4>) = —q. If one considers a body whose stationary 
configuration is specified by a single parameter, then the 
constraint M'((f>) = — q is enough to completely specify 
how the body changes, and therefore to completely de- 
termine its motion. This is the case in the weak-field 



1 Note that equations are to be applied on the near-zone 

field configuration, whereas equation JTJ refers to the background 
(far-zone) field configuration. In the paper these different con- 
figurations are distinguished notationally and their relationship 
is given precisely. 



2 It is sometimes stated that the baryon number is held fixed. 
This defines the sensitivity only for bodies with an identifiable 
baryon number that completely specifies the body. Examples 
where this definition does not work include black holes and finite- 
temperature stars. 
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and black hole cases worked out in section IV Bl below. 
For bodies specified by more than one parameter, the 
constraint M'(cj>) = — q only provides guidance, which 
must be supplemented by physical arguments or by find- 
ing slowly time-dependent solutions. 

The final component of the results provided by our for- 
malism is the equation for the perturbation to the smooth 
background spacetime produced by the small body. We 
find for a general theory of our class that the charges ap- 
pear in a universal way as delta-function sources for the 
linearized field equations, equations (|117p - (|119[) . below. 
Written out for Brans-Dicke theory, this becomes 



[V C V C 



(i) 



■'ab 



-8tt 



2oj 



(M + 2q<p)6 A (x,z(T))dT (5) 



8tt / MuaUbS^x, z)dr, 



(0) 



where superscript (1) denotes linearization, the sources 
on the right-hand side are constructed from the back- 
ground, and £ ao is the operator £ a t, = 0G o b+V c V c 0.g a f,— 
V a V b </> - STrT [ a f, with G ab the Einstein tensor and T* 1 
the scalar-field stress-energy, given in equation (|185[) . be- 
low. The integral is over the worldline of the body, 
denoted 7, and £(4(2;, 2/) is the "covariant delta func- 
tion" , given in a coordinate system by 8i(x, x') = S(x fl — 
X ' fl )/V~9- Once the mass ilf, charge q, and worldline 
7 have been determined, equations §5§ and ([5]) allow the 
associated waveform to be determined, making contact 
with far-field observables. These equations agree with 
linearized versions of those obtained by Eardley [6] 

Equations (HJ-© give a simple demonstration the na- 
ture of our results using the example of Jordan frame 
Brans-Dicke theory. While the new formulae @ and (O 
are certainly of interest in this classic theory, the main 
power of the present approach comes from its generality. 
Indeed, expressions for H}-© are given for an arbitrary 
theory of our class (equations (|115l) - (|123p 'l. and it is a 
simple matter to convert these expressions into useful 
forms for a given theory. We have performed this proce- 
dure not only for Jordan frame Brans-Dicke theory but 
also for generalized scalar-tensor theories in both Jordan 
and Einstein framesjfl for Einstein-Maxwell theory, and 
for the Will-Nordvedt vector-tensor theory. In most cases 
we have also evaluated the charges for some known so- 
lutions and worked out the corresponding predictions of 
the theory. 



3 Eardley's equations are only formally defined, since delta- 
function sources appear in non-linear equations. A natural regu- 
larization is to simply linearize his equations, and our formalism 
provides a rigorous justification for the resulting (well-defined) 
equations. 

4 One nice consequence of this analysis is the explicit demonstra- 
tion of how equation JTJ describes bodies in both frames, with 
the meaning of q and M changing in precisely the manner re- 
quired to ensure consistency. See equations l|157|l and l|158|l for 
the surface-integral formulae for the Einstein frame mass and 
charge. 



In section |TT] we review Lagrangian field theory and 
symplectic structure, establishing notational conven- 
tions. In section Hn] we review the framework and results 
of paper I and derive the new, general formulae for the 
mass and charges. In sections IIVI we apply the prescrip- 
tion to general relativity with a cosmological constant, 
while in sections [V] and IVII we consider specific scalar- 
tensor and vector-tensor theories. 

Our conventions are those of Wald Q ■ Early-alphabet 
lowercase latin characters a, b,... are abstract spacetime 
indices and greek characters are spacetime coordinate in- 
dices. When coordinates t, x l are used, mid-alphabet 
lowercase latin characters ... are spatial coordinate 
indices, while a denotes the time component. Wc 
also employ a spacetime multi-index notation using mid- 
alphabet capital latin characters, I, J,.... The meaning 
of these indices is described in detail at the beginning of 
section Hi A[ below. 



II. LAGRANGIAN FIELD THEORY 

In a series of papers spanning more than a decade 
Wald and coauthors developed a systematic and 
mathematically rigorous approach to Lagrangian and 
Hamiltonian field theory of a general n-dimensional co- 
variant theory, with applications to (e.g.) black hole ther- 
modynamics, conserved quantities, and stability. Here 
we review their basic approach and present the results 
that are relevant for our treatment of particle motion. In 
particular, we restrict their quite general results to our 
(still quite general) context of four spacetime dimensions 
and theories containing only scalar and vector fields in 
addition to the metric. The extension of our methods to 
higher dimensions and higher rank tensor fields is left to 
future work. 



A. Lagrangian and Symplectic Structure 

We consider a theory of a four-dimensional Lorentz 
metric g ab along with a collection of scalar and vector 
fields. However, for notational convenience, we will ex- 
plicitly consider only one scalar field </> and one vector 
field A a . (We indicate, where necessary, how each for- 
mula would be modified for the presence of additional 
fields.) We will denote the collection of all the dynamical 
fields {g, 4>, A} by ^0 We often employ a multi-index no- 
tation using mid-alphabet capital latin letters, I, J, K, .... 
Here a lowered multi-index represents the "natural" po- 
sition of the index on each field, <J// ~ {g a b, <j>, A a }, while 
a raised multi-index indicates the "conjugate" position, 
* / ~ {<? af \0i^a}- Free indices indicate the presence of 



5 Note that Iyer and Wald [131 use a lowercase i/> to denote the non- 
metric fields and a lowercase <b to denote the full set of fields. 
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an equation for each tensor field, while contracted indices 
indicate a sum over contractions of each tensor field, e.g., 
iffj^ff 1 — 9abg ah +'t> 2 +A a A a . When using multi-indices we 
do not distinguish notationally between coordinate and 
abstract indices. In the coordinate index viewpoint, ^ffi 
stands for "each component of each tensor field", with 
the potentially confusing aspect that the multi-index po- 
sition refers to the natural position, so that, e.g., the 
lowered-index iff / actually refers to the raised-index com- 
ponents of the vector field, A^. As seen below in equation 
(|10jl , we adopt the standard convention of varying with 
respect to fields in their conjugate index position (raised 
multi-index) , so that field equation operators are defined 
in natural position (lowered multi- index). Thus, we vary 
with respect to ^f 1 ~ {g ab , <fr, A a }, producing field equa- 
tion operators Ei ~ {E l £, E&\ E^ a }. We switch back 
and forth between multi-index notation and ordinary no- 
tation as necessary. 

A covariant Lagrangian field theory is specified by the 
choice of a scalar function L that is locally and covari- 
antly constructed @ from the dynamical fields iff and 
their derivatives, 

L = L[g ab ,cj ) ,A a ]=L^]. (7) 

Here and below, the use of square brackets indicates that 
the function considered is a local, covariant function of its 
arguments and their derivatives. (We make an exception 
for linearized fields, as described below equation (fl5)l .) 
In particular, we have the covariance property, 

a*L{ff} = L[ar**], (8) 

for diffeomorphisms a. In the passive viewpoint, this 
equation states that L transforms as a scalar, and we 
refer to L as the Lagrange scalar. 

The Lagrange scalar, as a scalar, cannot be integrated 
over the manifold without the specification of a volume 
element. In the Wald approach, one uses the natural 
volume element e a bcd associated with the metric g a b, us- 
ing a Lagrange 4-form Le a bcd as the fundamental object. 
The Lagrange 4-form may then be integrated over the 
(four-dimensional) manifold M without the need to spec- 
ify an additional volume element. However, for calcula- 
tions done in coordinates it is useful to introduce a fixed 
volume element e a bcd (also denoted d A x) and define \J—g 
to be the scalar function such that \/—ge a bcd — tabcd- 
(If e a bcd is taken to be the coordinate volume element 
associated with a coordinate system, then y/—g is given 
by the square root of minus the determinant of the met- 
ric components in that coordinate system.) In this case 
it is customary to work with a Lagrange density y/~gL, 
which is to be integrated with respect to the fixed volume 
element e a bcd — d A x. 

Many of the key results of the Wald approach to La- 
grangian field theory are most simply derived and ex- 
pressed in the language of differential forms |l3j]. How- 
ever, it can be notationally cumbersome to retain the 
indices associated with forms, and the results are not 



naturally expressed in a form that lends itself to calcula- 
tion in coordinates. While we will in some cases employ 
the differential forms approach in reviewing the results of 
the Wald approach to Lagrangian field theory, we will pri- 
marily use dualized versions of the quantities defined by 
Wald and collaborators (preferring scalars and vectors to 
4- forms and 3- forms). We will use the term "Lagrangian" 
to refer to the choice of L in a general way, while we will 
be careful to distinguish between the Lagrange scalar, 
Lagrange 4-form, and Lagrange density when necessary. 

Given a choice of Lagrangian, the action for the theory 
is written as (recall our notation e a bcd = d x), 

S = J L^g~ d A x. (9) 

The precise formulation of a variational principle based 
on this action requires a careful formulation of a region 
of spacetime and assumptions about the behavior of vari- 
ations on its boundary. Following jl 31 ] , we will avoid any 
associated delicate issues by defining all quantities of in- 
terest relative to the Lagrangian, without dealing with 
the action directly. 

We now introduce the notion of a variation of the dy- 
namical fields $ . For a smooth one-parameter family of 
field configurations ^(A) (not necessarily satisfying the 
equations of motion), we denote the background config- 
uration by VP (that is, *ff = \&(A = 0)) and the perturbed 
configuration by S^ff (that is, 8^ = (#a*(A))|a=o)- Tne 
variation of a function F(*ff) of the fields is then defined 
by considering the one-parameter family F{if!{\)). Iyer 
and Wald [l3] demonstrate that the variation of the La- 
grangian may always be written 

5{^/=gL) = + V o 0° [*,<$*]) , (10) 

where (as explained above) the square bracket nota- 
tion implies functions constructed locally and covariantly 
from the arguments and their derivatives. 

The tensor functions E] [iff] are called the field equation 
operators of the theory, and solutions iff to Ej [iff] = are 
called electrovacuum solutions. (For fields iff that do not 
satisfy the field equations, we would regard £7 ['J'] 7^ 
as a source; see discussion in section III C\ ) The vector 
function <d a is called the symplectic potential, defined up 
to an identically conserved current. Its antisymmetrized 
variation is called the symplectic current uj a , 

V^g w a [*, <5i*, S 2 ^f] = S 2 (V^g 6 a [*, - 1 -H- 2. 

(11) 

The symplectic current may be used to construct the 
Hamiltonian of the theory in the manner described in 
13]. However, for our purposes we require only the fol- 
lowing useful property, 

J^g~ V a Lu a = S 1 {y/=g Ei) <5 2 * 7 - 1 o 2. (12) 

This property follows from equations (|11[) and (fT0|) . to- 
gether with the fact that for a function F of the dynam- 
ical fields we have 8182F — 525\F by equality of mixed 
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partial derivatives. (Note that the calculation is consid- 
erably simpler if done in the language of forms, using 
Le a bcd along with the hodge duals of 8 a and u> a . In this 
case one may use the fact that the exterior derivative is 
independent of the metric and thus commutes with varia- 
tions.) Using S^J^g = -(l/2)g ab 8g ab y^g, equation (fl"2")) 
may equivalently be written as 

V a u a = -^g ab 8 1 g ab E I 8 2 ^ 1 + S 1 E I 8 2 y I - 1 o 2. (13) 

In particular, if the background configuration ^ is a so- 
lution of the field equations, then we have 



V a w Q [<J, 8^, 8 2 V] = E} 1 ' ft* ]<5 2 tf ' - EY ] [8^8^ 



(14) 



where we have defined the linearized field operators E\ 
in the standard way, 



(i) 



E { p[8V] = SEl 



(15) 



Note that Ey' is covariantly constructed from \1/ and 
S^f (and not 8^ alone). We make this one exception to 
our rule that the bracket notation indicates a function 
covariantly constructed from the displayed arguments. 
Here and below, a function notated with a superscript 
(1) will depend on a background field configuration ^ in 
addition to the displayed arguments. 

Equation (|14p has three important implications. First, 
it shows that when Si^ and S 2 ty satisfy the linearized 
field equations off of an electrovacuum solution ^, then 
one may construct a current uj a that is covariantly con- 
served on the background. Second, it shows that the 
linearized field operator off of a background solution ^> 
is self-adjoint (according to the definition given in [HI). 
Third, it shows how, when the linearized field operator is 
acting on one field and contracted with another, one may 
"move" the operator onto the second field at the cost of 
a total derivative. This type of operation occurs in the 
method of integration by parts, and will be of direct use 
to us in obtaining our surface-integral formulae for the 
mass and charges of a small body. 

In section IIIII and below we will assume that the field 
equation operators Ei are of second differential order. In 
this case, the linearized field equation operators E^p may 
be written 



Ef'lSV] = A/j[*]5* J + BJj[*]V 5* J 
+ Cf 7 [*]V a V h( 5* J , 



(16) 



for some A/j, BJ 7 , and Cfj covariantly constructed 
from the background configuration \& and its derivatives. 
Comparing with equation (fl4"|) and using the Leibniz rule, 
it follows that we may always choose the symplectic cur- 
rent to have the form 

uj a = T> a IJ [^}S 1 ^ I 8 2 ^ J 

+ F?j[V] (<Ji* J V 6 J 2 * J - V 6 rf<5 2 * J ) , (17) 



for some D° 7 = D} 7 and ~Ffj covariantly constructed 
from the background configuration and its derivatives^ 
Inherited from the corresponding freedom in O a , one al- 
ways has the freedom to modify the symplectic current by 
the addition of an identically conserved current. While 
we will assume in some intermediate calculations that uj a 
takes the form of equation (fTT|) . the final results (|120|) - 
(|123p are valid for any choice of ui a . (The invariance is 
explicitly ensured by Stokes' theorem.) Equation (fl7|) 
holds only for theories with second-order field equations. 
For higher derivative theories, an analogous form would 
be obtained with higher derivatives of the perturbations 
appearing on the right-hand side. 



B. Bianchi Identity 

Now consider the variation induced by a one-parameter 
family of diffcomorphisms (i.e., consider an infinitesimal 
diffeomorphism) . In this case the variations of the dy- 
namical fields are given by C^ 1 for a smooth vector 
field £. By the covariance property ([5]), the variation of 
the Lagrange scalar is simply 8L = C^L, and we have 



-g ab £ 6 g aD L + C^L 
=s(V„£ a L + £ a V„L) 



(18) 

(19) 
(20) 



However, from equation (|10j) we have 

8(y/=gL) = V=9 (E^C^ 1 + V a 6 a [*,£^]) , (21) 
and combining equations (PZUl) and (|2"Tj) gives 

V a (6 a [*, C^} - = -E^jC^ 1 . (22) 



Equation ([22]) shows that current 8" — £ a L is conserved 
if ^> satisfies the field equations or if £ generates a sym- 
metry of the dynamical fields (or both). In the langauge 
of Noether's theorem, this current is simply the Noether 
current associated with the local symmetrjlj of infinitesi- 
mal diffeomorphsims. However, we may obtain a current 
that is conserved even when EjC^ 1 7^ by noting that 
C^ 1 is linear in £ and its first derivative, so that the 
Leibniz rule may be used to convert this term into a term 
linear in £ as well as a total derivative. To perform this 
operation we work with the individual fields rather than 
the multi-index notation. (For the analogous calculation 



One may straightforwardly attain the relationship of these quan- 
tities to A/ j, BJj, and CJj by direct computation. Indeed, one 
may directly compute the symplectic current from the linearized 
equations in this manner. However, it is usually easier to proceed 
from the Lagrangian via equations I I1UI) and 1111 . 
By a local symmetry we mean a variation such that 8(^/—gL) = 
v / — g\7 a ct a for some ct a . Equation d 20 1) shows that ct a = £ a L for 
infinitesimal diffcomorphisms. 
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done in total generality, see reference [l7|-) In this case 
we have 



EjC^ 1 = E l a 9 b ] C i9 ab + eWjC s 4> + E^ a C^A a 
= £^ ] (-2V a £ b ) + E^C^a(t> 

+ E [A]a^b VbAa+AcVa ^ 

= V a [-2E^e + E^ a A c ^) 
+ ^ b (2V a E l a ^+E^V b (j) 

+ E^ a V b A a ~V a {EW a A b ) 
Combining equations (|22j) and (|25j) gives 



(23) 
(24) 



(25) 
(26) 

J£ = 6| - - E l a 9 b k b + E^ a A4 c (27) 

E^V b (f> + 2E^ a V [b A a] - V a E^ a A b . 

(28) 



with 



B a = 2V a E [ a 9 b ] 



where 0| reminds the reader that O a is evaluated on 
the perturbation C^ip 1 . Now suppose £ a (and hence J| ) 
has compact support, and consider the integral (using 
the volume element e = ^J—gcfix) of V a J a over a region 
containing the support region of J a . Stokes theorem im- 
plies that this integral must vanish for all £ a . Thus the 
integral of £, a B a vanishes for all £ a , and in fact B a must 
vanish. We thus derive that 



V a J£ = 
B a = 0. 



(29) 
(30) 



Equations (|29|) and (T30|) hold for arbitrary (not necessar- 
ily symmetry-generating) vector fields £ a and arbitrary 
(not necessary electrovacuum) field configurations \E ,/ . 
Thus we have an identically conserved current given 
by equation (f2"?| as well as a differential identity, 

V°4b - \^ a E^ a A b = E^ a V [a A b] - IeW V b <f>. (31) 

We refer to equation (|3"Tj) (equivalently equation (|30|) ) 
as the Bianchi identity for the theory. We emphasize 
that this equation holds independent of whether the field 
configuration {g, <j>, A} satisfies the field equations. If the 
theory contained additional scalar or vector fields, this 
identity would be modified by the addition of appropriate 
terms of the form present above. If the theory contained 
higher rank fields, there would be additional terms of 
similar structure; the explicit form of these terms is given 
in [l7|. 



8 Note that our J|* differs from that of Iyer and Wald [l3tl ■ who 
- ( a L. 



C. Non-Electrovacuum Solutions 

Before proceeding we make a few comments (and es- 
tablish notation) regarding solutions with matter sources, 
or "non-electrovacuum solutions" . We consider a La- 
grangian L with action S with field equation operators 
Ej, whose solutions Ej = are called electrovacuum 
solutions. Suppose, however, that one wishes to con- 
sider the case where additional matter fields are present. 
A common approach is to simply add an unspecified 
"matter action" Sm, so that the total action becomes 
S' = S + Sm ■ In this approach one does not specify ad- 
ditional dynamical matter fields to be varied with respect 
to, but rather leaves Sm as some unspecified function of 
unspecified extra fields and varies only with respect to 
the variables present in S. In this case the equations of 
motion become 



^ab ~ 

E m = 

E [A]a _ 



-1 SSm 

-1 SSm 

-1 SS M 
sf^g SA a 



1. 



;t„„ 



J a , 



(32) 
(33) 
(34) 



where SSm/S^ 1 is the functional derivative, which de- 
fines the matter stress-energy T ab , scalar charge density 
cr, and charge-current density J a with conventional nu- 
merical factors given above. As a point of interpreta- 
tion, note that the stress-energy of the fields W is not 
included as part of T ab , but rather appears in the met- 
ric field operator E 9 b (likewise for scalar charge and 
charge-current). For example, in Einstein-Maxwell the- 
ory we have E l a 9 b ] cx G ab - 8nT^ b M , where Tg M is the 
usual electromagnetic stress-energy tensor, so that T ab 
has the interpretation of the non-electromagnetic stress- 
energy. This is the origin of our name "electrovacuum 
solution" for solutions of the field equations with vanish- 
ing {T ab , a, J a }. 

Since not every classically admissible matter source 
will follow from a Lagrangian — and since the details of 
any matter model are irrelevant to our main results — we 
find it more natural to work at the level of equations of 
motion for the matter, rather than proceeding from a La- 
grangian. For our purposes, a "matter model" is a set 
of equations for some auxiliary matter variables (den- 
sity, pressure, fluid four- velocity, etc.) together with a 
prescription for constructing {T ab , a, J a } from these vari- 
ables. We then couple that matter model to our theory 
by solving 



E 



[g] 



1 



Tab 



ab 

EM = — (j 

fi[A]a _ ja 



(35) 
(36) 
(37) 



use J£ 



together with the equations for the auxiliary matter vari- 
ables. While we do not require our sources {T ab ,p, J a } 
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to follow from a Lagrangian, the reader more familiar 
with that approach may imagine that {T ab , p, J a } are ob- 
tained from a matter action from functional derivatives 
via equations (j3"2")) - (j3"4")) . Note that the Bianchi identity, 
equation (|31[) , implies that solutions will only exist when 



V Q T Qb - V a J a A b = 2J a V [a A b] + aV b 4>. (38) 



III. MOTION AND CHARGES 

We now review the methods and results of paper I, 
paying special attention to the appearance of the charges 
that characterize the small body and its motion. In the 
process, we will derive the new formulae for these charges. 
As in the derivation of the Bianchi identity, above, we re- 
strict for simplicity to a theory of a metric, scalar, and 
vector field. While most of the analysis of the section gen- 
eralizes immediately to a general theory, in some places 
the assumption of a single scalar and vector field plays 
an important role in the form of the expressions. In these 
cases we indicate how the derivation and results would 
change in the presence of additional fields. 



A. Framework and Assumptions 

The basic philosophy of paper I (derived from that 
of [IH) is as follows. We imagine a field configuration 
^> (the "physical configuration" ) which contains a finite- 
sized body that is small compared to the scale of variation 
of the external universe. (In particular, all scales asso- 
ciated with the body, such as mass, charge, etc., must 
be small compared to the scales of the external fields.) 
We then seek an approximate description of the physical 
configuration by attaching it to a one-parameter family 
^(A) to be Taylor expanded in A@ In order that the 
truncated Taylor expansion accurately approximate the 
physical configuration vp, we choose the family such that, 
as A — > 0, all scales associated with the body approach 
zero. Since Taylor's theorem guarantees an asymptotic 
approximation as A — > 0, our perturbation series in A 
should accurately approximate the physical configuration 
as long as the scales of the body are small compared to 
the scales of the external universe. 

Of course, to implement the above procedure one must 
identify suitable assumptions on the one-parameter fam- 
ily that reflect the satisfaction of the above physical crite- 
ria. Our approach, described in detail in jlLll9j. involves 



9 For smooth one-parameter families, this viewpoint is equivalent 
to the more heuristic approach of splitting the physical field con- 
figuration into a background and a perturbation and expanding 
in the perturbation. A treatment of small bodies will usually 
entail singularities on a worldlinc, in which case the equivalence 
is less evident and we find the formal mathematical approach 
significantly clearer. 



the requirement that certain "near-zone" and "far-zone" 
limits exist and are related to each other in a smooth way. 
The main output of these assumptions is an assumed 
form of the field perturbations in coordinates t, x % that 
have the interpretation of being centered on the body 
(i.e., the worldline x % — 0, called 7, describes the lowest- 
order position of the body). We also introduce the dis- 
tance coordinate r — yJb~ijX % xi as well as n l — x l /r, de- 
noted by n when representing a direction on the sphere. 
We use superscript (m) to denote an mth order pertur- 
bation, 



*i(A) 
* 7 (A) 



A* 



(i) 



X^ 1 + 0(X 2 ), 



(39) 
(40) 



where, as is evident from these equations, we do not use 
the unperturbed metric to raise and lower indices of the 
perturbations. This convention harmonizes more natu- 
rally with the variational viewpoint, in terms of which 
the symplectic current — and hence our effective mass and 
charges — are naturally expressed. However, it will be 
useful to have a map that relates vfrt 1 ^ to and we 

define qu and q IJ by 



1 



(41) 
(42) 



In matrix language equation (|4ip reads 



o (1) 

9qb 



1 c 




1 



(l)o'b' 




while equation (14"2"T) reads 



(43) 





_ g {a)aa'g(0)bb' 


A (0)b' § a' 

(44) 

The "components" of qu and its inverse q IJ may be read 
off of the matrices in equations (14"3"|) and (j4"4")l . respec- 
tively. Note that this raising and lowering convention 
does not allow for the raising and lowering of a single ab- 
stract index on the metric perturbation , but only for 
the raising of both indices at once. We will use this rais- 
ing and lowering convention for the majority of the paper. 
However, when displaying expressions for the symplectic 
current, we will considerably shorten the length by rais- 
ing and lowering individual indices with the background 
metric (and abandoning the convention described above). 
In order to avoid confusion, in all such cases we will use 
a special symbol for the perturbation, such as /i M „ in- 
stead of g$ , and in all such cases we will explicitly note 
the convention. Note also that, as evidenced in equa- 
tions ([4"U|) - (|4"4")) . we do not bother to maintain a constant 
index position with respect to perturbative superscripts 
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(or other non- index modifiers); we happily write both 
*< 0) and <j(°) z . 

At finite A, we impose the field equations only for 
r > XR (for some R), leaving the interior of this world- 
tube arbitrary; this ensures that our results hold inde- 
pendent of the internal composition of the body. The 
perturbation series, then, is only defined for r > 0. How- 
ever, we assume that the background fields are smooth 
(the body "smoothly disappears"), so that we may in- 
clude the curve r = in the background field configu- 
ration. We mark the value taken at r — with a hat, 
writing 



H> I (t,x i )=4> I (t) + 0(r). 



(45) 



Here and below the 0(r n ) symbols are at fixed t,n. We 
will always set the metric on the worldline 7 to rj^ by 
coordinate choice, so that equation (|4"5)l reads 



(46) 
(47) 
(48) 



9® = »W + 0(r) 
<£(°) = 4>(t) + 0(r) 
A*® 1 * = i"(t) + 0{r). 



Our assumptions imply that the first perturbations take 
the form 



(1)/ _ 



ai(t,n) 
r 

a 1 (t, ft) 



+ 0{r°) 
+ 0(r°) 



(r > 0) 
(r > 0), 



(49) 
(50) 



i.e., the perturbations diverge as 1/r near 7. We have 
introduced the field a/ to characterize the time and an- 
gular dependence of the divergence. As is evident for the 
above equations, we raise and lower indices on ai with 
the map qu rather than the unperturbed metric. For 
the individual fields, we denote a/ = {ajfj,, a^ M }, so 
that we have 



= 


a{fi(t, n) 
r 


- 0(r°) 


(r > 0) 


(51) 


0« = 


aM(t,n) { 
r 


- 0{r°) 


(r > 0) 


(52) 




a^(t,ft) 


+ 0(r°) 


(r > 0). 


(53) 



Equations (|40[) - (|53|) describe the far-zone background 
and perturbations $W. In this far-zone limit, the 
body shrinks to zero size/mass/charge and disappears 
entirely from the background field configuration, leaving 
behind the smooth "external universe" 'J' ' . However, 
the body leaves behind the preferred worldline x % — as 
the place where it "disappeared to" , and the field per- 
turbations become singular on this worldline. 

We also consider an alternative, near-zone limit, which 
is designed to keep the body at fixed size. While the far- 
zone limit represents the viewpoint of an observer far 
from the body, where its effects are small, the near-zone 



limit represents the viewpoint of an observer in the vicin- 
ity of the body, where the effects of the external uni- 
verse are small. To define this limit we first introduce 
rescaled coordinates x l — X~ 1 x t and (for some time in) 
i = A _1 (i — to). An observer at fixed scaled coordinate 
moves correspondingly closer to the body as it shrinks in 
size, with all measured scales reducing to zero. However, 
if we "change units" by rescaling the fields, 



9fiv — A g^v 

4> = <t> 

= XA*, 



(54) 
(55) 
(56) 



then the scales associated with the body (which shrank 
to zero in the far-zone limit ]now remain fixed, and a 
finite A — > limit is achieved^ Note that the powers of 
A in the definitions of the rescaled fields precisely cancel 
the powers of A that arise from a coordinate transfor- 
mation to scaled coordinates. Thus the components of 
the scaled version of a field can always be computed by 
simply plugging in for the new coordinates, 

g^{X;t ;i,x l ) = g^(X;t = t + Xi,x l = Xx l ) (57) 
0(A; t ; t, x l ) = 4>(X; t = t + Xi, x l = Xx l ) (58) 
A>*{X;to;i,x l ) = A" (A; t = t + Xi, x l = Xx l ). (59) 

We denote scaled coordinate components by placing a 
bar on the coordinate index, so that the above equations 
relate components in scaled coordinates of scaled fields to 
corresponding components of the original fields in origi- 
nal coordinates. 

Having introduced the scaled coordinates and fields, 
the scaled limit is simply the A — > limit, 



mi 



lim g^(X]t Q \t,x l ) 



<f> (U) (t ;t,x l ) = lim (j)(X;t ;t,x z ) 
A^ttolt,^) = lim Af"{X;t Q ;i,x l ), 

A— >0 



(60) 
(61) 
(62) 



where the limit is taken at fixed scaled coordinate. Since 
we work only at r > A-R (for some R), the limiting fields 
, J'( ) are defined only for f > R, i.e., only for sufficiently 
large f. Thus, while we apply the term "near-zone" to 
this limit (and refer to as the "near-zone background 
field configuration"), it should be borne in mind that 
these fields correspond to a near-zone description only 
of the body's exterior. (Indeed, for this reason these 
fields were referred to as the "body exterior fields" in 
1].) Determining the relevant near-zone fields in a given 
application requires matching to an interior solution, as 
discussed below equation (|6"5)l . below. 



10 Note that defining rescaled fields with conjugate index position 
would require different powers of A. For simplicity we will not 
define such fields in this paper, working only with scaled fields 
in the "natural" index position of equations H54I I- H56I1 . 
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An immediate consequence of equations (j57l) - (j59|) is 
that the near-zone background fields are stationary (in- 
dependent of t), and, coupled with our more general as- 
sumptions, a further consequence is that the fields admit 
the large-f expansions, 

#(*o;^) = ^ + ^^ + o(i) (63) 

^)(t„;^) = 0(to) + ^^ + o(l) (64) 

AMfa;* 1 ) = >(*o) + a[A] "f^ + o (1) , ( 6 5) 

Note the dual role played by the hatted quantities as the 
value of the far-zone background metric as the position 
of the particle and the asymptotic value of the near-zone 
field configuration. Similarly, note the dual role played by 
the al*l as the 1 fr divergence of the far-zone perturbation 
at time t and the 1/f correction to the asymptotic value 
of the stationary near-zone field configuration computed 
at time t = to. 

1. The Near-zone Background Field Equations and a 
Subtlety for Non-Scale-Invanant Theories 

An important point regarding the near-zone back- 
ground fields is that that these fields do not in gen- 
eral satisfy the full field equations, but rather obey those 
of a "scale-invariant sub-theory" . For the scaled variables 
at finite A, we define 

^[^^^^^[A^^A- 1 ^] (66) 

for each field This defines an operator on field config- 
urations as a function of A. If the equations scale 
homogeneously, [X 2 g, 0, \~ 1 A] = X N E^ [g, <f>, A] for 
some integer TV, then the operator agrees (up to pow- 
ers of A) with the original operator E^\ and the scaled 
fields satisfy the same equations as the original fields. If 
not, then in general the scaled fields will satisfy a differ- 
ent equation. In particular the relevant operator for the 
background scaled fields is 

£ [ * It0) &,Ml = Jim X N E^[X 2 g,^X^A], (67) 

A->0 

where the integer N is chosen in each theory in order 
to ensure that the limit exists and is non-trivial. Finite- 
A satisfaction of the field equations (for r > A-R) then 
implies 

EW(%<-°\fi a ),AW]=Q, (f>K). (68) 

We will refer to i^K ) as the near-zone exterior field 
operator. The near-zone fields obey this equation for 
sufficiently large f, while for f < R the fields are not 
constrained by the formalism (indeed, we have not even 
defined "I^ ) in this region). Instead, the formalism makes 



reference only to properties of the exterior fields 
(through the charges associated with the 1/f falloff at 
large f). In this way we obtain results independent of 
the detailed composition of the body (depending only on 
the charges). However, in order to apply these results 
one must determine the exterior fields and charges asso- 
ciated with a given body type, which necessarily involves 
a continuation of into the interior f < RVH 

How should this continuation be done? The details 
will depend on the particular theory under considera- 
tion, but a basic guide is the physical interpretation of 
the large-f region as representing the "buffer zone," far 
enough from the body that its fields may be approxi- 
mated as monopole (1/f) fields, but still close enough 
that the external universe's fields may be approximated 
as constants. In the physical (finite A) solution, equation 
(|68p will be approximately satisfied in this zone, and in 
particular the approximation will break down as one ap- 
proaches the body. To mimic this behavior in our per- 
turbative solution, one should model the body interior 
and surrounding region using the full, non-electrovacuum 
field equations and, at some distance away, attempt to 
interpolate to a solution of the restricted equations 
that has the asymptotics of equations (|6"3"|) - (|6"5|) F^1 This 
ensures that the properties of the near-zone fields for 
large f accurately reflect the buffer zone properties of 
the body one desires to model. 

The reason for the necessity of this awkward inter- 
polation is that our mathematical description maps the 
physically finite buffer region to the infinite coordinate 
region of sufficiently large f . This allows the charges to 
be precisely defined in terms of asymptotic f —¥ oo be- 
havior, but at the cost of the satisfaction of the full field 
equations, so that the interpolation described above is 
required. An alternative, physically equivalent approach 
would be to consider solutions to the full field equations, 
but evaluate the surface integral formulae only in a buffer 
region. The reader that prefers this viewpoint may re- 
gard all of our surface integral formulae as holding for 
approximate spheres in an approximate buffer zone. 

If the theory is invariant under the scalings of equation 
(|6"6")) . then these subtleties do not arise. To determine the 
near-zone fields of a body one simply considers the full 
non-electrovacuum field equations coupled to the matter 
model of choice, demanding stationarity and the asymp- 
totic boundary conditions given by equations (|M|) -(|55 |) . 



We remind the reader that R does not represent the radius of the 
body, but is simply an arbitrary number designed to ensure that 
we impose the vacuum field equations only outside the body. 
Note that the full field equations may not even have solutions 
satisfying these asymptotically flat boundary conditions, indi- 
cating the general necessity of doing this kind of transition for 
modeling an isolated body. 



10 



2. Summary of Consequences of Assumptions 

Returning to the multi-index notation, we may com- 
pactly summarize the important consequences of our as- 
sumptions as 



*j. 0) =*/(t) + 0(r) 



(i) _ ai(t,n) 



(0) 



*j(*o) 



+ 0(r°) 
ai{t ,n) 



(69) 
(70) 

(71) 



where satisfies the (electrovacuum) field equations, 
satisfies the linearized (electrovacuum) field equa- 
tions (off of vfr(°)) for r > 0, and satisfies the near- 
zone exterior field equation, equation (|68l) , for sufficiently 
large f. (Each ^K ) is defined individually in (|54"1) - (|6"2"1) .) 
We will always set the "hatted metric" to r\ by coordi- 
nate choice, i.e., = {77, cj), A^}. We remind the reader 
that indices of the perturbation (|70p are not raised and 
lowered with the background metric g(°\ and that the 
near-zone background fields are always considered with 
natural index position , rather than conjugate posi- 
tion 



B. Point Particle Description 

Our above assumptions proceed entirely from consider- 
ations regarding finite, extended bodies, as spelled out in 
detail in 111, 

El- 

However, the resulting far-zone pertur- 
bations are singular at r — (going like 1 jr) , which sug- 
gests a connection — within linearized theory — to effective 
point particle sources. (Readers familiar with the distri- 
butional equality V 2 (l/r) = — 4ir5(x) may readily attest 
to this fact.) This connection was exploited in 0, [l9| . 
leading to the result that the far-zone perturbations are 
sourced by effective point particle sources characterized 
by a finite number of charges. We now repeat/improve 
the derivation of this result, while going beyond the treat- 
ment previously given in that we provide a formula for 
the effective charges. 

The far-zone perturbations, equations (j5Tj) - (j53|) , sat- 
isfy the linearized field equations for r > 0, 

4 5 i (1) l5 (1 \0 (1 U (1) ]=O (r>0) (72) 
BW(%W,^),iW]=0 (r>0) (73) 
ijMW/^W 0« A«] =0 (r>0), (74) 

which in the multi-index notation reads simply 

£jW[*( 1 )]=0 (r>0). (75) 

Since the perturbations diverge only as 1/r, however, 
they are locally integrable, and hence may be straight- 
forwardly promoted to distributions defined everywhere 
(including r = 0). 



To define distributions we take our space of test func- 
tions to be C°° compact support functions and per- 
form integrals using the natural volume element e = 
yf—gd^x associated with the background metric . 
(We write this element as 



-g, even though \J — g(°) 
might be a more accurate notation.) We define distribu- 
tions with both raised and lowered multi-index. Raised- 
index distributions act on lowered-index test functions 
fi = {fjfv ; /'^, /'^l, while lowered-index distributions 
act on raised-index test functions. We raise and lower 
indices on distributions by raising and lowering the test 
function using the map gu introduced in equations (j4Tj) - 
(|44[) . For example, given a lowered-index distribution T>i, 
the raised-index version is defined as 



(V I J I ) = (V I ,q IJ fj), 



(76) 



where the notation (a, b) indicates the action of a distri- 
bution a on a test function b. This convention guarantees 
consistency with the convention of raising and lowering 
ordinary functions with qrj, and allows us to blur the 
distinction between raised and lowered index distribu- 
tions. Thus we will often write (V, /) instead of (Vj, f 1 ) 
or(^,//). 

Returning to the task at hand, we associate the dis- 
tribution ^a 1 ' with the locally integrable function 
by 



~gd A x 



(77) 



Since the linearized field equations are linear, their oper- 
ation on a distribution is automatically defined by 

(£«[§«],/) = J J—gd±x E^im^ 1 , (78) 

where the self-adjointness of E^p has been used (see 

remarks below equation (TUT) ). Since Ej 1 ^ acts on the 
smooth function /, the integral is well-defined, and we 
may equivalently write 



(£W [§«],/) = lim 



^dtxEPlfi&W*. (79) 



However, since the integral now ranges only over points 
where r > 0, we may make use of the fact that the per- 
turbations satisfy the linearized field equations on 
that domain. To do so, we use equation (fl4|) to "move 
the field operator from the test function to the field", 
giving 



[$(1)], /) = lim 



-gd x 



(4 1) [* (1) ]/ I + V B w [*W ) /,*< 1 



(80) 



The first term vanishes by equation (|75j) . while the second 
may be converted to a surface integral by Stokes theorem, 

(£« [#«],/) = lim / \/^hd 3 x(N a LJ a [^ \f,^ i y 

(81) 
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where \/—hd a x is the 3-volume element associated with 
the induced metric h on the surface r = e (a timelike 
world tube for sufficiently small e), and N a is the nor- 
mal vector. To proceed further we adopt our locally 
Minkowskian coordinate system, equation (|46l) , where we 
have g^ v = r)^ v + 0(r). In particular we have ^/—hd 3 x = 
r 2 (l + 0(r))dfldt, N a = O(r), and Ni = -n t + 0(r). Fur- 
thermore, from the form of u a (equation (|17|l ) together 
with the fact that and / are smooth, we have 



= / 



•0(l/r). 



(82) 



Thus equation (|8ip becomes 



(%f) = - lim / AM n^[* (0) ,/,* (1) ] 

r— >0 J 

(83) 

--fW'WU&f**^' (84) 
= - | dt(/ J F&)| 7 | dfi ^ - n j( A , (85) 

where in the last line we have used the form of the per- 
turbations, equation (1701) . However, since the action of 
i^W [VfrC 1 )] depends only on the test function evaluated at 
r = (i.e., on 7), the distribution E^[^/^] is in fact a 
delta function, 



with 



Nj = 47rF} J 7 [*(°)] 



da \ 



(86) 



(87) 



where the double angle brackets represent an average over 
n. Once the symplectic current (and hence Fjj[^]) is 
known, equation (|57|) gives an explicit formula for the co- 
efficients of the delta function effective source in terms of 
the background solution Vl^ ) and the coefficients a 1 (t, ft) 
of the 1/r divergence of the perturbations "I^ 1 ). Since 
these coefficients may be determined in the near-zone 
(see equation (|7ip). equation (|57|) contains a prescription 
for determining the charges of a body from its fields. In 
practice, however, it is easier to work directly with the 
symplectic current uj a rather than the function Fjj that 



was defined in equation (|17l) . In light of equation 
and the antisymmetry of uj a we may rewrite equation 
(157)1 as 



Nj(t) = lim / r 2 dfl 

r^Oj 



riiUj 



a 1 (t, n) 



cept the Jth, where it equals one[[3 Note that the in- 
dex position is important here; it is required that the 
raised- index components f 1 be assigned these values (i.e., 
Si 1 ^ 1 = S 1 j 7^ Si^i), corresponding to the "conjugate" 

positions of the field values, f 1 = , f l \}. 



Equation (|88|) gives a convenient prescription for deter- 
mining the charges from the background (far-zone) so- 
lution ^j -* and the coefficients a 1 , which may be deter- 
mined from the near-zone field configuration, equation 
(f7l"j) . corresponding to the body of interest. Note that 
in many cases of interest the symplectic current will de- 
pend only on ^ = ^W| 7 (i.e., there is no dependence 
on derivatives of the background configuration), which 
also appears in the near-zone field configuration. In this 
case equation (|88[) may always be rewritten as a formula 
making reference only to the near-zone solution 



Application of Bianchi Indentity and Derivation 
of Formulae for Motion and Charges 



Equation (JHH) shows that the first-order perturbations 
have effective delta-function source. We may gain sig- 
nificantly new information about this source through the 
use of the Bianchi identity for the theory. To proceed 
further we restrict explicitly to a theory of a single scalar 
and vector field (in addition to the metric), although it 
is straightforward to continue in total generality (see sec- 
tion V of paper I). After linearization, the Bianchi iden- 
tity (f3"Tj) becomes 



V°£^ 1(1) - ~V £^(i)«4°) 



(89) 



where V is the derivative operator associated with the 
background metric . This identity is linear in the 
perturbations ty^ 1 ' and therefore holds automatically on 
the distributional versions In particular, applying 

the identity to equation (|86p indicates that solutions will 
exist only when 

V*(JVWtf(a:*)) - l -V^{N^5{x l ))A^ 

(90) 

where we have dropped the superscript (3) on the delta 
function for readability. Working out the consequences 
of this equation is simplified by use of Fermi normal co- 
ordinates (e.g., [2(|) about the worldline x l = 0. In these 



(88) 



S 1 j indicates that, in deter- 



where the notation f 1 
mining the Jth component of TV, one should choose the 
test function f 1 which equals zero in all components ex- 



13 Although a good test function must be of compact support, the 
integral in equation (1881 ranges over only a finite domain and 
one may instead work with the function f 1 = 5 j. 
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coordinates the metric takes the form 

9$ = V„u - tptuOttyx* + 0{r 2 ), (91) 

where tf, — (—1,0,0,0) and a l (i) are the spatial com- 
ponents of the four-acceleration of the worldline, a M = 
(0,a l ). The Christoffel symbols in the Fermi coordinate 
system are given by 



r 



oo 



1 i0 



0(r), 



(92) 



with all other symbols vanishing on the worldline (i.e., 
T% = 0(r) except as indicated above). In terms of the 
Fermi coordinates equation (|90[) becomes 

= - 2d N$ + 2a* N$ - d N^°A 

- aiN^Ao + N^do<t> - 2N a d [a A 0] } 



2N„ 



[a] 



N\ A] A 



■} 



(93) 



= 5(x){ - 2d N$ + 2a j N\f + 2oiN$ - d^N^At 

- aj N^ j Ai + N^di<j> - 2N^ a d [a A i] } 

+ 8 : >5(x){2Nlf-N lA \A i y (94) 

where the superscript (0) on the background fields 
and A^ has been dropped for readability. To determine 
the implications of these relations we integrate against a 
test function and use the fact that its value and derivative 
are arbitrary. This shows that the terms proportional to 
djS(x) must separately vanish, while the terms propor- 
tional to S(x) must equal the spatial divergence of those 
proportional to djS(x). From the vanishing of the djS(x) 
terms we have 



2 N W=nI A %, 



(95) 

which provides constraints on the charges iVj. In par- 
ticular, the anti-symmetric part of equation (|95[) shows 
that N^ A \ must be proportional to A[ \ and the symmet- 
ric part then determines in terms of the constant of 
proportionality. We thus have for some %(t) that 



Nl A] = X A 



N, 



[g] 



^xAiA„, 



(96) 
(97) 



where we remind the reader that superscript (0)'s have 
been dropped. Thus the Bianchi identity in fact reduces 
the number of independent charges from fifteen [Nj = 

{N$,NW», ATI*!}) to just four, which may be taken to 

be Njjjj, N^°, x, and N^. We give these charges the 
suggestive names, 



Arte] — I m 
iV oo — 2 

= -q 
N [A]0 = e 

X = e. 



(98) 

(99) 
(100) 
(101) 



The remaining relations in equations (|93[) and (l94l) (those 
proportional to the value, rather than derivative, of the 
test function after integration) then give 

d M = -A d Q e - ^edo(A i A i ) - qd <t> (102) 

Ma t = d Q eA t ~ -e(A j Aj) + d (eAiA ) 

+ 2ed [0 A l] + qd^ (103) 

where we again drop the superscript (0) on the back- 
ground vector field A^ ^ and scalar <^°) for readability. 
Equations (fT02]) and (|T03l) give the mass evolution and 
equation of motion for the small body in Fermi normal 
coordinate form. We can also rewrite equation (1861) in 
light of equations (IMl) - (|10ip . giving 



oo 

[<?](!) 
iO 

[g](i) 
E l<t>]W 

E [A](1)0 
E [A](l)i 



-M5(x) 
\eA^Af6{x) 

\eA^Af5(x) 

-qS(x) 

e5(x) 

eA {0]t 8{x) 



Finally, the formulae for the charges ([55)1 become 



-M = lim 

2 r-y0 



-q = lim 

r^0 



lim 

r-s-0 



eA l = lim 

r->0 



r 2 d£l UiUj 1 
r 2 dfl niLO 1 
r 2 dil 7ijo/ 
r 2 d£l UiUj 1 



(0) 



-Wo, 0,0} 



* (0) , -,{0,1,0} 
r 

* (0) ,- l {o I o,0 



(0) 



-, {o,o,5V} 



(104) 
(105) 

(106) 

(107) 
(108) 
(109) 

(110) 
(111) 
(112) 
(113) 



where it is understood that e vanishes if A 1 vanishes. 
(Since only the product eA l appears in our equations of 
motion, we may as well make this simplification.) In the 
case that A 1 does not vanish, we can provide an explicit 
expression for e by dotting with Aj and using the form 
(|82f, giving 



lim 

r->0 



r 2 d$l niUj 1 



(0) 



i 1 r S\Ai 
-,{0,0,-^— 
r 1 A3 A, 



(114) 



We remind the reader that in equations (|110|) - ()114j) the 
index position is important for the third argument of w°; 
the symplectic current is viewed as a function of raised- 
multi-index variations Si^ 1 and 62^ , and one must send 
62^* (rather than 62^1) to the prescribed values above. 
For example, to compute the mass, one sends 829^ = 
5Vo> corresponding to S 2 g^ = -6°S°. 
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D. Covariant Presentation of Results 



Equations (|102[) - (| 114[) constitute our results for gen- 
eral theories, given in Fermi coordinate form. We now 
rewrite these results in covariant form, in terms of the 
four-velocity u a of the worldline and the projection or- 
thogonal, P ab = g ab + u a u b . We have for the motion 
that 



M 



e (u»A b ) 



u b \/ b u a = P ab \u c V c eA b 
\eV b (P cd A c A d ) + u c V c (eu d A d A b! 



(115) 



+ 2eu c V [c A b] +gV b< />}, 
for the mass evolution that 

u a V a M = u a ^-u b V b eA a + ieV a (P bc A b A c ) - gV^j , 

(116) 

for the linearized fields that 

E mi) = - I q 5 4 (x,z(T))di 

E [A](l)a 



(117) 



(u a [e + e(A c u c )} + eA a ) 5 4 {x, z)dr (118) 



E 



[am 



(u a u b [M - e(A c u c ) 2 ] + eA a A b ) 



x 64 (x, z)dr, 
and for the charges that 



M = 2 lim / r 2 dn muz' 

r->0 



q = — lim / r 2 dQ riiUj 1 

r->0 



$ ( ",-,{uV,0,0} 



* (0) , -,{0,1,0} 



e = lim / r d^t riiio 1 

r->0 



e = lim / r 2 d£l riiU) 1 

r->0 



*<o) « ,{o,0,-^} 



• {0,0, A^} 



(119) 

(120) 
(121) 
(122) 
(123) 



where Aj^ = P ab A b / (P cd A c Ad) is a vector pointing along 
the spatial projection of A a . (It is understood that if 
spatial projection vanishes, P ab A h ~ 0, then the charge 
e vanishes and equation (|123l) docs not apply.) In all 
of these equations we have dropped superscript (0) on 
background fields. While equations (I115p - (|119[) are fully 
covariant, we have left equations (|120|) - (jl21[) in a form 

valid only in local inertial coordinates, gffi — r/^ + 0(r). 
(Note in particular that — (1,0,0,0) and n % = x l /r.) 
Since any computation of the charges will likely use such 
coordinates, we see little reason to covariantize further. 
We remind the reader that index position is important 
here; see discussion below equation (|114l) . For example, 
in computing the mass one must send <52<7 MI/ = u^u^, 
corresponding to 629^ = —u^u^. 



We now indicate how the results (|115p - (|123l) gener- 
alize to theories containing additional fields besides the 
single scalar and vector field considered here. First, it 
was shown in paper I that analogs of equations (|115[) - 
(|119P exist for a theory with arbitrary types of tensor 
fields, though the explicit form was not obtained. As 
for the formulae for the charges, new to this paper, our 
general formula (|88[) remains valid for a theory with ar- 
bitrary types of tensor fields, where the multi-index runs 
over each component of each field. However, the appli- 
cation of the Bianchi identity to simplify this formula, 
leading to a smaller number of charges (just four in the 
scalar- vector-tensor case), must be redone to determine 
analogs of (|120l) - (|123[) for a general theory. Thus, for 
an arbitrary theory, the status is that analogs of our re- 
sults (|115l) - (|123p are known to exist, but would take a 
complicated form that has not yet been written down. 

An exception to this statement, however, is the case 
where the theory conditions additional scalar fields only. 
Following through the calculations of section lTlTl for a the- 
ory of a metric, a vector field, and a collection of scalar 
fields 4>a, it is easy to establish that equations (|1 15|) - (|123[) 
generalize as follows. First, there is a new scalar charge 
qA associated with each field 4>a- To generalize equations 
(I115P and (|116|) , one simply adds an index A to each inci- 
dence of the scalar field 4> and its charge q, turning these 
terms into sums over the collection of fields. For exam- 
ple, the last term of equation (|115p . qP ab \7 b (f), becomes 
qAP ah y b 4>A, where there is a sum over A. Equation (|1 17[) 
is modified by creating an additional copy for each field: 
one adds an index A to the field <f>, charge q, and field 
equation operator £#](!). Equations (|TTg|) and ([TTO1) are 
unmodified. Finally, there is a copy of equation (|12ip for 
each new charge, where the second argument of the sym- 
plectic current contains a non-zero entry only for the slot 
associated with the field whose charge is under consider- 
ation. Equations (|120p . (|122p . and (I123[) are unmodified, 
except in the trivial sense that the symplectic current has 
additional arguments for the additional scalar fields, into 
which zeros are to be inserted. 



E. Usage of Results 

We now summarize how the results of this section may 
be applied in a specific theory. First one must obtain the 
near-zone fields for the body type of interest as a func- 
tion of internal parameters (e.g., baryon number, cen- 
tral pressure, etc.) and the local external field values, 
as described in section fill A II One may use any asymp- 
totically Minkowskian coordinate system (t,x l ). Next 
extract the asymptotic 1/r fields in that coordinate sys- 
tem, i.e., determine the a/r terms defined by equations 
([63| - (p5)l . Then plug the asymptotic fields into equations 
(|T20 | - (fT23l) with = (1, 0, 0, 0) to determine the charges 
of the body as a function of the internal parameters and 
external field values. (See also equations (|110p - (|113l) . 
where coordinate components are written out explicitly.) 
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After selecting a background spacetime, plug the charge 
functions into equations (|115[) and (|116l) , where the value 
of the background configuration {g, A, 4>} on the world- 
line is used as the external field value for the charges. For 
bodies with a single internal parameter, equation (|1 16[) 
will fix the evolution of this parameter, and the coupled 
system may be solved for the motion. For bodies with 
more internal parameters, evolution laws for the param- 
eters must be determined by physical arguments or by 
finding near-zone solutions with slowly varying asymp- 
totic fields. However, once evolution laws (consistent 
with (|116p ) for the paramters are fixed, one may solve 
equation (|115l) to determine the motion. The motion 
and charge evolution may then be plugged in to equa- 
tions (|117j) - (|119|) to determine the perturbation caused 
by the body. Note that it may be convenient to elimi- 
nate the step of explicitly extracting the a/r near-zone 
fields by rewriting the surface integrals in the near-zone, 
as done for each specific theory treated below. 



IV. GENERAL RELATIVITY WITH A 
COSMOLOGICAL CONSTANT 

We now work out the consequences of the general pre- 
scription for mass, charge and motion in the example of 
general relativity with a cosmological constant. The La- 



grangian is 

L = ±-(R-2A)., 

which gives field equation operator 

M 



16nEW[g]=G ab +Ag ab 



(124) 



(125) 



(126) 



and symplectic potential 

lGirQ a [g,Sg] = 2g b M a Sg b]c . 

The symplectic current is then given by 

327Tu a [g, h, j] = j ab V b h + j be V a h bc - 2j bc W b h c a 

- 3 V a /i + jV b h ab -h+*j. (127) 

Here we have set h ab = Sig ab and j ab = 5 2 g ab , and we 
raise and lower indices on these quantities with the back- 
ground index, so that h ab = -8\g ab and j ab = -d 2 g ab - 
Since there are no fields besides the metric, the only non- 
vanishing charge is the mass M. To determine a formula 
for this charge, equation (|120[) (equivalently (|110|l ) in- 
structs us to consider the symplectic current with h^ v — > 



aT/r and f v -> 



superscript [g] on ajfi, since our theory has no additional 
fields. Note also that /i M „ = — a Mt/ , since we do not use 
the background metric to raise and lower indices of a^.) 
Computing u/ and keeping only the most singular (1/r 2 ) 
terms gives 



(We have dropped the 



327rcji 



9, ' 



di( 



O-hk 

\ r 



O 



(128) 



The mass is then given by equation (|120l) to be 

(129) 

Taking into account the appearance of a^„ in the near- 
zone background (equation (I63|l ). we may equivalently 
write 

M = iib & / f2dn nl ~ ' (130) 

which is recognized as the famous ADM mass formula Q ■ 
Thus, the charge M is nothing but the ADM mass of the 
near-zone background spacetime. It is notable that the 
ADM formula for M emerged, rather than some other 
formula (such as that of Komar [2l|) that is equivalent 
given the stationary field equations. Since the ADM ap- 
proach uses the Hamiltonian formulation, which itself is 
intimately connected to symplectic structure (e.g., [HI), 
one might suspect a deep and general connection between 
that approach and ours. However, we have found no 
such connection, and it may be pure coincidence that 
the ADM formula emerges. 

To determine the mass M of a given body using equa- 
tion (|130[) . one must determine the near-zone background 
metric. This procedure was discussed at the conclusion 
of section UlIAl The first step is to identify the near- zone 
exterior field operator. Using the scale-invariance of the 
Einstein tensor, at finite A we have 



= G ab [g] + Ag ab = G ab [g] + X 2 Ag ab 



(131) 



Letting A — > in scaled coordinates, we then obtain 

G pp [5 (0) ]=0 (larger), (132) 

which is satisfied for f larger than some R (see equation 
(l68l) and discussion below). Thus the near-zone back- 
ground metric is a solution of the vacuum Einstein equa- 
tions in the large-f region. Equation (p3"j) indicates that 
we should seek stationary, asymptotically flat solutions; 
and, following the discussion of section lTlI A 11 these solu- 
tions should be matched to interior solutions relevant for 
the particular matter model under consideration. The in- 
terior solutions should satisfy the full non-electrovacuum 
field equations, 



G^[.9 (0) ]+A5^ = 8^r pp 



(interior) (133) 



where T ab is the stress-energy tensor of the matter fields. 

Thus, after choosing a matter model one should at- 
tempt to find a solution of equation (1133[) that interpo- 
lates in a reasonably smooth way to a stationary, asymp- 
totically flat solution of (|132j) . This provides the near- 
zone background metric for the body under considera- 
tion, whose mass M may then be computed by equa- 
tion (|130l) . (This is just the ADM mass of the solution.) 
Note, however, that while in principle the cosmological 
constant term in equation (I133P should be included, in 



15 



practice its contribution should be negligibly small. This 
is because the requirement of a region surrounding the 
body where the solutions match approximately implies 
that the body must be small compared to the scale set 
by the cosmological constant, in which case its effects 
should be negligible everywhere in the near-zone. Thus 
one may simply consider the field equations of ordinary 
general relativity (with no cosmological constant) when 
determining the mass of a particular body type. In more 
general theories, however, this type of argument may not 
be possible, and one must in principle use the full field 
equations near the body. 

Having interpreted the charge M, we may now make 
use of the results for the motion and metric perturba- 
tion, equations ()115[) - (|121j) . In particular, the motion is 
geodesic, 

u b V b u a = 0, (134) 

the mass is conserved, 

u a V a M = 0, (135) 

and the linearized perturbation h is sourced by a point 
particle of mass M, 

G ab\h] + A/i ab = — J u a u b 8 A (x, z)dr. (136) 

To summarize, the prescription for general relativity 
with a cosmological constant is as follows. First, specify 
what type of body is being considered by giving a matter 
model and determining the relevant stationary, asymp- 
totically fiat solution to equations (|133[) and (|132[) (the 
cosmological constant term in (|133|) may be neglected). 
Compute the ADM mass of this spacetime and assign 
this value to the parameter M. Given a background so- 
lution of the vacuum Einstein equation with cosmological 
constant, the motion of that body will be geodesic, and 
the perturbations produced will satisfy equation (I136p . 

V. SCALAR-TENSOR THEORIES 

For a theory of the metric and a scalar field, (|116[1 - 
(I115P become 

u a V a M = -qu a V a <j) {0) (137) 
Mu b V b u a = q (V 0)afc + u a u b \ V 6 (o) . (138) 

The main new subtlety in the scalar-tensor case (com- 
pared to general relativity) is that the charges q and M 
are not conserved. Indeed, the mass varies according to 
equation p37[) . while the evolution of q is simply not 
constrained. The lack of deterministic evolution is a re- 
sult of imposing field equations only outside of a world- 
tube: While in general relativity satisfaction of the vac- 
uum field equations in this region suffices to determine 
the motion of a body, in a scalar-tensor theory we only 



learn that the motion is given in terms of two charges, 
one of which is not constrained at all by the theory. (See 
paper I for more discussion of this point.) Of course, if 
the region containing the body is "reinserted" at finite 
A for a given scalar-tensor theory and a given matter 
model, then one should have deterministic evolution at 
finite A and hence deterministic evolution within our per- 
turbation series. The precise deterministic evolution will 
depend in detail on the precise matter model considered, 
but our results indicate that this freedom will show up 
only in terms of a single charge q (the mass M being 
given by equation (| 13T[) once known initially). Thus, to 
complete the prescription of the motion of a body in a 
scalar-tensor theory, one must determine the value and 
evolution of q implied by specific matter models for spe- 
cific types of bodies. 

We now discuss this process. Each scalar-tensor theory 
will have a different set of interior and exterior field equa- 
tions, analogous to equations (|133l) and (|132j) . that the 
near-zone field configuration {g^°\ (j)^} should satisfy, 
with a smooth interpolation in an intermediate region. 
Comparing equations (|4l)|) - (|4"Tl) and we see that 

the f — > oo boundary conditions are that the metric must 
approach the Minkowski metric and the scalar field must 
approach the value of the far-zone scalar field evaluated 
at the current position of the particle. That is, for a given 
time to we seek stationary (t-independent) solutions of 
the near-zone equations such that 

^(io^-^V (f->oo) (139) 

(o) (io;x 4 ) $(t ) = ^\ t=to ^ =0 (f ^ oo). (140) 

How does one determine the appropriate such solution? 
The simplest method is to specify a matter model to- 
gether with physical assumptions (such as a fixed total 
baryon number) such that the solution will be uniquely 
determined for each possible external field cf). In this case 
the field configuration {g^°\ <f>' ^} becomes a determinis- 
tic function of the local external field value <j>. The mass 
and charge are then also determined uniquely as func- 
tions of 4> by equations (|157p and (|158p . and we write 

M = M(4>) (141) 
q = q{4>). (142) 

In this case we have u a V a M = M'(0(°))u a V Q (o) and 
equation (|137[) implies that 

M'{4>) = -q{$). (143) 

Equation (|143p gives a constraint that must be satisfied 
by any valid prescription for obtaining the charges {M, q} 
as a function of the external field value <j>. In light of the 
constraint we may interpret the charge q as represent- 
ing the sensitivity of the mass of the body to adiabatic 
changes in the local external 0-field value. The notion 
of a sensitivity was first introduced by Eardley [|| in 
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the context of a phenomenological approach to motion in 
Jordan frame Brans-Dicke theory, based on a point par- 
ticle action. Will and coauthors have subsequently ap- 
plied and expanded the approach, as summarized in |22j | . 
Damour and Esposito-Farese [23] adapted the Eardlcy 
method to a more general class of scalar-tensor theories, 
and furthermore related q and M to the asymptotic 1/r 
behavior of the fields. Our first-principles treatment has 
confirmed the appearance of a sensitivity when one con- 
siders bodies whose exterior field configuration is a strict 
function of the local external field value (a basic assump- 
tion of previous approaches, relaxed by ours). The sen- 
sitivity is normally defined by s = —d\ogM/dlog(j), and 
so the relationship to our q is given by 



M 



(144) 



We emphasize, however, that the relations (I14ip - (|142[) 
(and hence (|143[1 ) are not fundamental to scalar-tensor 
theory but rather are consequences of considering a body 
whose exterior fields are "state functions" of the of the 
local external field value (p. While this is surely a reason- 
able assumption for many bodies of interest, there is no 
principle that forbids more general bodies, whose exter- 
nal fields may be "path functions" of (j) or may change 
independently of </> due to internal processes. An exam- 
ple of the latter would be a star that cools or otherwise 
changes its internal structure on a timescale relevant for 
motion. (Since the stellar structure and evolution equa- 
tions may be drastically modified in alternative theories, 
this possibility is not obviously excluded for our Uni- 
verse, as described by alternative theories.) In this case 
one would instead have a prescription for determining the 
exterior fields {g^°\ 4>^} as functions of both </> and some 
prescribed internal dependence related to (say) a slow 
change in equation of state. That is, these fields would 
become explicit functions of the time to along the world- 
line, in addition to any dependence on external fields. In 
this way the mass and charge would become functions of 
both to and 4>, and one must respect the more general 
constraint (|137p , which becomes 



A. Einstein Frame Scalar-Tensor Theories 

We now consider general relativity with the addition 
of a Klein-Gordon field with potential V(<p), 



L 



R 

16^ 



1 



ff ab V a ^V b 0- V(<t>). 



(148) 



It is well known that a large class of Lagrangians can be 
written in this form by a field redefinition, a procedure 
normally referred to as "going to the Einstein frame" . In 
doing this procedure one normally assumes that the mat- 
ter was coupled to the original field variables in a stan- 
dard way ("minimally coupled"), so that the coupling 
to the new field variables becomes non-standard. This 
assumption amounts simply to a choice of non-standard 
matter model for the Lagrangian (| 148|) . so that many 
scalar-tensor theories are subsumed under the single La- 
grangian. In this section we derive the surface-integral 
formulae for the charges for the Lagrangian (|148[) and 
evaluate them for specific types of matter coupling. 
The field equation operators for the Lagrangian (|148[) 



Gab _ 

2 ab 



fJ 



16tT 

ab v Q v, 



<P-V'(cf>), 



(149) 
(150) 



where 



T ab = V a 0V h </> - -g ab g ccl V c 4>V d <i> + V(^)g ab (151) 

is normally referred to as the stress-energy of the scalar 
field. The symplectic potential is 



6 a = & a GR - V Q ^ 



(152) 



where Qqr is the symplectic potential for general rela- 
tivity, given by equation (|126|) . The symplectic current 
is then given by 



u a [{g, </>} , {h, a} , {j, (3}} - u% R \g, h, j] = 



P ( V a a + - W a (j 



h ah V b cj> 



a <H> P 



, (153) 



■£-M(to,$) = -q(t Q ,4>)u a V a 
dt 



For completeness we reproduce equations (|117|) - (|119j) 
for the body fields in the scalar-tensor case, 



E l<P)( 



E 



ab 



q 5 4 (x,z{T))dT (146) 



Mu a Ub&4,{x, z)dr. (147) 



where we view u as a function of k 

,,ab 



ab 



Sig ab , ] ah 



(145) $29 , ot = Si4> and P = <5 2 0, and use the background 
metric g ab to raise and lower indices on these quantities. 
(The last term on the RHS stands for an additional copy 
of all previous terms on the RHS, except modified by 
exchanging h with j and a with p.) The formulae for 
the charge M, equation (|120j) . instructs us let h^ v — > 
aWff/r, j» v -> u^u v = 8%8 v o, a -> /r and p -> 0. 
Keeping only the most singular (1/r 2 ) terms, we have 



>] 



J GR 



(154) 
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where oJqr is evaluated with the same arguments as uj 1 
(given explicitly by equation (|128p ). Thus the formula for 
the mass M is identical to that of general relativity (i.e., 
to the ADM mass of the near-zone spacetime). For the 
scalar charge q, equation (|121[) instructs us to let h^ v — > 
a^ u /r, f v -> 0, a -> a^/r and /3 -> 1. Keeping only 
the most singular (1/r 2 ) terms, we have 

>] / fl W\ / 1 



Thus the charge g is given by 



lim / r 2 dVLn l d, 

r->0 



(155) 



— ■ (156) 



Taking into account the appearance of affl in the near- 
zone background field, equation (164)) . we may rewrite this 
formula in terms of the near-zone background field. We 
display the result along with the (ADM) formula for M, 



M = —— lim 

lD7T r->oo 



— lim 



r 2 d£l n 



f 2 dn n% 



1 fag, 



(0) 
hi 



d i9kk 



(157) 
(158) 



To determine M and q for a particular body we must 
determine the near-zone fields g^ and <j>(°K Using the 
scaling properties of the field operators, it is easy to see 
that we have 

G^[5 (0) ]=87r(v^(°)Vp^°) 

- ^V 0) ^V a ^>V^°>) (159) 
g(o)&^_ V( g^ 0) = 0, (large f) (160) 

where V a is the derivative operator associated with g^ . 
These equations agree with the full field equations with 
the potential V{6) set to zero. On the other hand, the in- 
terior satisfies the full non-electrovacuum field equations 
(including the potential V((f>)), 



(161) 



V 



) = -a, (interior) (162) 



where T a b and a are the stress-energy and scalar charge 
density of the matter under consideration. To deter- 
mine the charges q and M one should find solutions of 
equations (|159p - (|162[) in the manner described in section 
IIII A 11 and plug these solutions into equations (|157l) and 
(|158|) . The motion is then given by equation (|138[) . and 
the far-zone body fields are given by (|146|) - (1147p . 



Ordinary Coupling: Einstein- Scalar Theory 



Consider the Lagrangian (|148|) with no potential, 
V(6) — 0, and in the case that "ordinary matter" is 



endowed with an intrinsic scalar charge density a in ad- 
dition to its stress-energy T^. (In a Lagrangian view- 
point, these would be defined by equations (|3"2"|) and (|3"3")l .) 
We refer to this theory as Einstein-Scalar Theory. We 
now investigate the meaning of the charges q and M in 
this theory. We therefore seek solutions to the near-zone 
equations (|159p - (|162p , which, given V{4>) — 0, reduce 
to equations (|161[) - (|162p where the sources T a \, and a 
have compact support. These solutions should be sta- 
tionary and approach {77, <fi} at large distances (see equa- 
tions ipg p -ipD )i ). 

We first consider a "Newtonian body" in the sense that 
1) only linear deviations from {77, <p} are allowed and 2) 
Toi and Tij are negligible. Following the standard ap- 
proach to the Newtonian limit in general relativity, we 
find that a gauge may be chosen where equations (I16ip - 
(EUg) become 



V 2 $ = 47175 
-a, 



V 2 a 



(163) 
(164) 



with p = Too and V 2 = S^didj, and with the fields given 
in terms of the potentials $ and a by 



9^ = f)iiv - 2$^ 
(f> = (f> + a. 



(165) 
(166) 



From the multipolc expansion we know that the solutions 
to equations (|163p - (|164p that are vanishing at infinity are 
given by 



pd*x )- + 



1 



(167) 
(168) 



Plugging equations (|165l) - (|168p into equations (I157P and 
(|158p then gives the charges for a Newtonian body to be 



M = / pd 3 x 
q = I crd 3 x. 



(169) 
(170) 



We thus see that for weak-field bodies the mass and 
charge are given by volume integrals of the mass den- 
sity and scalar charge density, respectively. Note that 
these densities are not conserved, but are related to each 
other by equation (|38p . This equation does not fix their 
evolution in time, which will instead depend in detail on 
the matter model adopted. Thus the Einstein-scalar the- 
ory displays a wide variety of motions of bodies, even in 
the weak-field limit. Indeed, all we can conclude is that 
the general equation of motion, equations (|137p - (ll38p . is 
obeyed for some (in general non-constant) charge q. 

It is also straightforward to consider black holes. The 
well-known "no scalar hair" theorem [24]-[27j implies that 
the scalar field must be constant, 6 — 6. and hence the 
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metric must be that of the Kerr family. In this case the 
mass M agrees with that of the Kerr solution and the 
charge q vanishes. Thus the mass is conserved and the 
black hole moves on a geodesic. 



2. Conformal Coupling with no Intrinsic Scalar Charge: 
Generalized Brans-Dicke Theory 

In Brans-Dicke theory and its cousins, one assumes 
that matter is minimally coupled to a ( "Jordan frame" ) 
metric that is conformally related to the ("Einstein 
frame") metric of our Lagrangian (|148|) . This is rep- 
resented by considering a matter action of the form 



S 



M 



d A x n- 4 (ct>)V=9L M [i}M,n 2 (<t>)9 



ab] 



(171) 



where ipM is some matter field. Notice that the matter 
action depends on (f> only through the conformal metric 
^ -2 <7ah, which is identified with ordinary matter by the 
relationship 



^ord 

1 ab — 



-2 



SS 



V-det(ft- 2 (0) 5 ^) 5{W(<p)g-b) 
-2ft 4 (0) SS m 



~g 5(n 2 (<f>)g» b y 



(172) 



(173) 



where T° d represents the stress-energy tensor of "ordi- 
nary matter", i.e., the stress tensor one would normally 
assign to a matter model within general relativity. (This 
is the "Jordan frame stress tensor".) The stress tensor 
and scalar charge density of our Einstein frame theory 
are then related by 

-2 SS M 



Sg ab 
-1 SS M 
f=g S(f> 



« 2 T Q r 



Q- 3 n'g ab T° r b d , 



(174) 
(175) 



and so the non-electrovacuum equations are simply 

G ab ~ SnTlf = 8^ 2 T° b rd (176) 
.g ah V a V b = n- 3 n'g ab T£ d . (177) 



For a weak-field body (see equations (|163[) - (|164l) and text 
above), we then have 

(178) 
(179) 



v 2 $ = 47rrr 4 p 
v 2 a = -n~ 5 n' P , 



where S7 and its ^-derivative S7' are evaluated at the back- 
ground value (f> = (f>. Here p = S7 T 00 is the ordinary mass 
density, i.e., the time-time component of T ord coordinates 
Q,x^ where the metric Sl~ 2 g a b equals the Minkowski met- 
ric. Following the procedure of the previous subsection, 
the charges are given by 



M = SI 4 / p d 3 x = n^Mt 



q = si 5 o' P d 3 x = sr 2 si'M f 



(180) 
(181) 



where Mg — J P d 3 (flx) is the total ("Jordan frame") 
baryonic mass, i.e., the spatial integral of the time-time 
component of T ord in coordinates where the where the 
metric fl~ 2 g a b equals the Minkowski metric. From equa- 
tions (|180[) and (|18ip we see that the charge-to-mass ra- 
tio of Newtonian bodies depends only on the asymptotic 
scalar field value, 



(182) 



so that the motion of a Newtonian body is independent 
of its composition. (Note that this charge-to-mass ra- 
tio is precisely that required for geodesic motion in the 
metric Q,~ 2 g a b, consistent with the fact that matter is 
by assumption minimally coupled to that metric.) The 
constraint dM/d(f> = — q shows that the baryonic mass 
Mb is independent of the asymptotic scalar field value, 
in agreement with what is found in the Jordan frame in 
the following section. 

For black holes in Brans-Dicke type theories, the con- 
clusion is the same as in the simple Einstein-scalar the- 
ory: the no-hair theorems imply that that the metric is 
Kerr and the scalar field is constant, so that the scalar 
charge q vanishes, the mass M is constant in time, and 
the motion is geodesic. For bodies "intermediate" be- 
tween weak-field and black holes, such as compact ob- 
jects like neutron stars, one must solve the field equations 
and determine the mass and charge via the surface inte- 
grals (I157[) - (1158[) . One interesting example is provided by 
reference [28|, which identified a class of Einstein- frame 
scalar-tensor theories without potential in which compact 
objects can develop large scalar charges even in the case 
of small self-gravity. Their notion of scalar charge, inher- 
ited from the general work of 23] , agrees with our notion 
up to normalization. Accounting for a field redefinition 
relating their scalar field to ours, the precise relationship 
between our charge q and the charge w of 28] is given by 
q = \/4tt(jj. 

While we have performed all analysis for general con- 
formal factor, we note for completeness that in the orig- 
inal Brans-Dicke theory the conformal factor is given by 

Q 2 ((f)) = e 75=T3. 



B. Jordan Frame Brans-Dicke Theory 

Although the Brans-Dicke theory can be analyzed in 
the Einstein frame with no loss of generality, for com- 
pleteness we present the Jordan frame version. This ex- 
ample also illustrates nicely the generality of the force 
law, equation (|138|) . since it applies in both frames, with 
the meaning of q and M differing in precisely the way 
that ensures consistency of physical predictions. The La- 
grangian for Jordan frame Brans-Dicke theory is 



ujg 



ab V Q V^ 



(183) 
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where a; is a constant. The associated field equation op- 
erators are 



l&ir^g = 0G ab + U^g ab - V a V fc - 8tt7* 
IQttE^ =R- uj(j)- 2 {\7(t>) 2 + 2w0~ 1 D0, 



(184) 



with 



8nT, 



If = ^ (v a ^V b - X - (V0) 2 . 9a ^ , (185) 



where (Vcj)) 2 = g ab V a <fiV b (j> and U<j> = g ab \7 a V b <f). Note 
that these equations may be combined to yield 



□0 = 



16tt 
3 + 2w 



4>E 



(186) 



It is usually assumed that matter does not couple fun- 
damentally to the scalar field, E^ = —a = 0. We then 
have for the non-electrovacuum equations that 



ct>G ab + U4>g ab - V a V b = 8n(T™ + T ab ) 



n</> 



16tt 



g -Lab 



3 + 2oj- 

The symplectic potential for the theory is 

167r8 a = 2cj)g bc W la 5g b]c + W b <j>Sg bc 
- V a 4>g bc 5g bc -2uj4>- 1 V a ' 

and the symplectic current is 



(187) 
(188) 



(189) 



Taking into account the appearance of afi and in 
the near-zone background fields, equations (1531 and (j6"4"|) . 
we may rewrite this formula in terms of the near-zone 
background field configuration, 



M = — |- lim / r 2 dn n 

lD7T r->oo 



^ M? - a,® 



-2a 



(193) 



For the scalar charge q, equation (|121[) instructs us to let 
W a^/r, 0, a a^/r and /3 -> 1. Keeping 

only the most singular (1/r 2 ) terms, we have 



fc'[{^},{^-},{0,l}] 



r / \ r 




(194) 



Thus the charge g is given by 
1 



q = lim / r d£ln l 

I6n r->0 



+ 5i 



>]' 

'00 



2u}(f>- 1 d i i^ — 

J9]\ 



*fe*s 



(195) 



16^K[{ 5 , 0}, {ft, a}, {j, /?}] - ^) = --hf b V b cf> 

+ jV a a + u}(j}- 1 /3(-V a <j)h + 2h ab V b <j) + 2V a a) 
+ /3(W a h - V b h ab ) - {h, a} O {j, (3}, (190) 



Taking into account the appearance of a^l and in 
the near-zone background fields, equations (j63]) and (|64l) . 
we may rewrite this formula in terms of the near-zone 
background field configuration, 



where, as before, we view a; as a function of h ab — Sig ab , 
j ab = 52g ab , a = Si4> and /3 = S2(f>, and use the back- 
ground metric g ab to raise and lower indices on h and 
j. For the mass M (equation (|110[) or (|120|l ) we let 
W a ^ v /r, ->■ u»u v = 5^5% a -> a^/r and 
P — > 0. Keeping only the most singular (1/r 2 ) terms, we 
have 

lSnw'ttg, </>},{— - — }, 0}] = 




oil 

r 



(191) 



The mass M is then given by equation (| 120[) as 

i /■ T f 

M = lim / r 2 dn n l 

16n r->0 




(192) 



q = lim / r dfln 1 



(0) 



, Q ^(0) Q -(0) 



(196) 



In the Einstein frame, the mass formula involved only the 
metric, while the scalar charge formula involved only the 
scalar field. In the Jordan frame, both formulae involve 
both fields. 



1. Weak-field Bodies 

To put the formulae for the mass and charge to use we 
must find stationary solutions to the non-electrovacuum 
equations, (|187[) - (|188p . (Since the Brans-Dicke theory 
is invariant under our scalings, we do not need alterna- 
tive equations for the near-zone exterior.) We begin with 
weak field bodies. Assuming only linear deviations from 
a background {rj^, 4>} (where is a constant) and taking 
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Tqq = p and T K ; = 0, we find (following the original work 
[J]) that a gauge may be chosen where 



2 -87175 



3 + 2cj 



v^2 o l-i 5 + 4a; 

3 + 2w 



V 2 7y = 



(197) 
(198) 
(199) 



where the fields are given by 

g^u = V + Jfiv - -V^V^lap (200) 
= + a (201) 
The solution satisfying our boundary conditions is 



2 Mb , Q f I 



> x 2 + a; Mb ^ n ( 1 

.9oo = -1 + 40 — — +0 

3 + 2ui r 



(202) 
(203) 



9(j = 4u ( 1+4 -.i^^) + o(-L), po*, 

where we have defined the total baryonic mass Mb, 

M b = J pd 3 x. (205) 

Plugging this solution (as g(°\ 4>^) into equations (|193|) 
and (|196|) gives the mass and charge as 



M = M B 
9 = 0. 



(206) 
(207) 



This reproduces the fact that weak-field bodies in Brans- 
Dicke theory behave identically to their counterparts in 
general relativity: the mass is just the baryonic mass, 
and the charge vanishes, ensuring geodesic motion and 
constancy of M = Mb- We remind the reader that q 
is related to the sensitivity s by equation ()144|) . so that 
the sensitivity vanishes for weak-field bodies. Since the 
sensitivity is dimensionless and vanishes when the mass 
equals the baryonic mass Mb, it is often referred to as 
the gravitational binding energy per unit mass. 



2. Black Holes 

For a black hole in Brans-Dicke theory, the metric is 
given by the Kerr metric and the scalar field is equal to its 
asymptotic value 0. Adopting asymptotically isotropic 
coordinates, we have 



5oo 



9%j 



-i + ^ + ofl 

r \r z 



2M, 



GR 



o 



(208) 

(209) 
(210) 



where Mqr is the mass parameter of the black hole, equal 
to the ADM mass of the metric. Evaluating the mass 
and charge (equation (|193[) and (|196|) ) of the black hole 
solution gives 



M = <j}M GR 
1 

q=--M GR 



(211) 
(212) 



The sensitivity of the black hole, equation (|144[) , is given 
by 



a = 1/2, 



(213) 



in agreement with previous work. From the constraint 
dM/d<f) = — q (equation (11431) ). it follows that the mass 
parameter Mqr evolves with the external scalar field 
value by 



Mqr = Mq0~ 1//2 , 



(214) 



where Mo is some constant independent of 0, which 
might be viewed as the "intrinsic mass" of the black hole. 
In the weak-field limit, 1 ' 2 has the interpretation of be- 
ing proportional to the local value of Newton's constant 
G. Equation (I214p says that it is the combination GMqr, 
rather than the mass parameter Mgr, that remains con- 
stant for a black hole. Combining equations (|211[) and 
(I214P gives 



M = M^ 



(215) 



Thus the mass of the black hole scales as G -1 / 2 . If one 
constructs a "modified Planck mass" using G rather than 
G, then the black hole mass scales as the Planck mass. 
This scaling is sometimes assumed a priori, leading to an 
independent argument that the sensitivity of a black hole 
is 1/2. 



C. Generalized Jordan Frame Brans-Dicke Theory 

A simple generalization of Brans-Dicke theory studied 
recently by [29| is given by the following Lagrangian, 



L = ^U>R-^)g —— 



M(0) 



(216) 



Here the Brans-Dicke parameter uj has been promoted 
to a function of 0, and one has added a potential term 
denoted M(0). While these changes affect the equations 
of motion of the theory, they do not alter the symplectic 
potential or the symplectic current. Thus, the formulae 
for the mass and charge are unchanged, and equations 
(11571) and (|158l) hold for this theory as well. However, for 
a non-zero potential M(0), the theory will not be scale- 
invariant and the construction of the near-zone metric 
will involve the subtleties discussed in section UlI A 11 In 
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particular, to construct the near-zone metric one must 
interpolate between a full solution in the interior and 
a solution of the equation without the potential term in 
the asymptotic region. (Alternatively, a full solution may 
be used, with equations (|157[) and (| 158[) evaluated at an 
intermediate radius.) More details of this procedure are 
given in section UlI A II 

VI. VECTOR-TENSOR THEORIES 

As examples of vector-tensor theories we will con- 
sider Einstein-Maxwell theory and the theory of Will and 
Nordvedt [3(| ■ This latter theory offers a simple example 
of a non-Maxwell vector-tensor theory with second-order 
field equations. Many other such theories can be con- 
structed [Hj]; however, the simple Will- Nordvedt theory 
suffices to illustrate the features of our approach and the 
general properties of motion in the vector-tensor case. It 
is straightforward to repeat these calculations for a more 
general vector-tensor theory and derive the analogs of 
equations ([2TS1)- ([2"M1) , below. 

Vector-tensor theories whose vector field is constrained 
to be unit timelike are known as Einstein- Aether theories 
[3TM33I ] . Since a Lagrange multiplier scalar field is used 
to enforce the unit constraint, these theories are of a 
scalar- vector-tensor type according to our classification. 
However, while the Lagrangian is perfectly covariant, the 
scalar field does not behave appropriately in the buffer 
zon^]3 and must be eliminated as a dynamical field before 
the methods of this paper would apply. This elimination 
does not fit naturally into the framework developed in 
sections ILTllIIIl and a discussion of Einstein- Aether theory 
will be saved for a future publication. 

A. Einstein-Maxwell Theory 

The Lagrangian is 

L=-±-(R~g ab g bd F ab F cd ) 1 (217) 

1D7T v ' 

where F ab = 2V[ Q ^4fc]. The field equation operators are 

WttE^ = G ab - 2 (F a c F bc - \g ab g bd F ab F cd ) (218) 
E [A]a = ^_v b F bQ (219) 

The symplectic potential is 

© Q = ® a GR ~ ^-F ab 5A b , (220) 



For example, for a static linearized solution with a point mass 
source in the Einstein- Aether theory with C2 = C4 = 0, ci = 
— C3 (in the notation of |33|). the Lagrange multiplier field itself 
is a delta function, instead of having the desired 1/r behavior 
possessed by the other fields. Our formalism assumes that all 
fields in the Lagrangian have this behavior. 



and the symplectic current is 

uo a -u a GR = A (h ac F cd /3 b + F a d h bd fi b + 2/3 b VI a a b l) 

-{h,a}^{j,j3}, (221) 

where h ab = S l9 ab , j ab = 5 2 g ab , a a = 6 } A a , /3 a = 6 2 A ai 
and indices on these quantities are raised and lowered 
with the background metric. The Einstein-Maxwell the- 
ory has an additional gauge symmetry that may be used 
to simplify the calculations. This gauge symmetry gives 
rise to an additional charge-conservation identity for the 
theory, V a E^ a = 0. If this identity is used at the lin- 
earized level during the calculations of section UlI C[ one 
finds that e = and it a V a e = 0, i.e., that the hatted 
charge vanishes and the ordinary charge is conserved. 
(These calculations are done explicitly in paper I.) To 
provide formulae for the non- vanishing charges we use the 
symplectic current in equations (|110p - (|113[) and rewrite 
in the near-zone, following the pattern established with 
previous calculations for scalar-tensor theories. A short 
calculation then yields 

M = m / f2dn nl - (222) 

e=— lim / f 2 dfl n%A^\ (223) 

47T f-KX> J ° ^ ' 

We thus see that the mass M is simply the ADM mass, 
while the charge e is simply the electric charge as defined 
by Gauss' law. Using that e = and u a V a e = 0, the 
equations of motion (|115j) - (|116p reduce to 

Mu b V b u a = 2eu b V [b A a] (224) 
u a V a M = 0. (225) 

That is, the body follows the Lorentz force law and the 
mass (as well as the charge) is conserved. We thus re- 
cover the standard result for the motion of bodies, while 
connecting M and e directly to near-zone solutions via 
the surface- integral formulae (|222|) and (|223j) . We also 
recover the standard point particle charge-current and 
stress-energy, as may be seen from equations (|118[) - (|119p 
with e = 0. 



B. Will-Nordtvedt Theory 

The Lagrangian for the Will- Nordvedt theory [H |30| 

is 

L=-^(R + V a A b V a A b ). (226) 

The field equation operators are 

16irE [A]a = -2V b V b A a (227) 
IGkE® = G a b - 1 1 g ab V c A d V c A d 

+ V c {A c V (a A b) - A [o y c A b) - A (a V b) A c ) 
+ \7 a A c \7 b A c + \7 c A a \7 c A b , (228) 
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and the symplectic potential is given by 

167r(9 a - Q GR ) = Sg bc {A c V a A b + A c V b A a - A a V b A c 



28A b V a A b . 



(229) 



Taking an antisymmetrized variation to obtain the sym- 
plectic current, one obtains a long expression. Since only 
terms in which a derivative operator acts on a variation 
can affect the formula for the mass and charges, we only 
display these terms. The formula is 

167r(w a - u a GR ) = ~2n ab (5 b - 2n^ Jbc A c + A a JbcK bc 
-{h,a}^{j,p} 

+ (terms with no V on {h, a, j, /3}), 

(230) 

where 

Kab = V a a 6 + \A c {\J a h bc + W b h ac - W c h ab ), (231) 

and, as in the Einstein-Maxwell case, we have set h ab — 
Si9 ab , 3 ab = S 2 g ab , a a = 8 1 A a and (3 a = S 2 A a , and 
indices on these quantities are raised and lowered with 
the background metric. Furthermore, indices on K a b are 
raised and lowered with the background metric. (We note 
that K ab is just the variation of the derivative of the field, 
K a b = Si(V a A b ).) Following the pattern of calculations 
established previously in the paper, we find near-zone 
formulae for the mass, charge, and hatted charge are, 

M = — lim / r 2 dVL n % (d k g k i - dig kk 

107T r-voo J \ 

+ 2{A°V A t - A.VaA - A°V 4 A°)) (232) 



e = — lim / r 2 dn n,V l A° 



-1 

8ir r- 

e = — lim I r 2 dfl n.A.V'A 1 . 

8?T r^oo A k A k 1 



(233) 
(234) 



Here we have dropped bars and superscript (0) for con- 
venience; however, as in analogous formulae elsewhere, 
one should use the near-zone background configuration 
{g < f b \A^ a } in these formulae. (In particular V is the 
derivative operator associated with g(°\ Note also that 
the near-zone configuration is stationary; covariant time- 
derivatives in these formulae could be eliminated in terms 
of partial spatial derivatives of the metric.) It is under- 
stood that when A 1 = then the charge e vanishes. (Re- 
call that A a gives the asymptotic value of the vector field, 
lim^oo A^ a = A a .) Equations (P2>(|2"M)) provide for- 
mulae for the charges associated with a body in the Will- 
Nordvedt theory. Once the field equations have been 
solved for a given body type, these equations, combined 
with equation (|115p with q = 0, may be used to determine 
the motion of that type of body. Since the hatted charge 
is in general non- vanishing and the ordinary charge is not 
conserved, the Will-Nordvedt theory — and presumably 
other more general vector-tensor theories — allow a much 
richer phenomenology of motion than Einstein-Maxwell 
theory. This phenomenology remains to be explored. 
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